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Abstract 

The small sub-barrier tunneling probability of nuclear processes can be dramatically enhanced by 
collision with incident charged particles. Semiclassical methods of theory of complex trajectories 
have been applied to nuclear tunneling, and conditions for the effects have been obtained. We 
demonstrate the enhancement of alpha particle decay by incident proton with energy of about 
0.25 MeV. We show that the general features of this process are common for other sub-barrier 
nuclear processes and can be applied to nuclear fission. 
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Tunneling in nuclear processes has been a subject of study for many years since this is a 
substantial mechanism of nuclear decay and nuclear reactions, including nuclear fission and 
fusion (see, for example, and references therein). The recent interest in understand- 

ing the processes of under-barrier tunneling has been stimulated by the 

calculation of bremsstrahlung radiation in alpha decay when the alpha particle is moving 
under the barrier j^. In this letter we consider another feature of quantum tunneling: the 
possible enhancement of nuclear decay due to interactions with low energy charged parti- 
cles. This enhancement has the same origin as the tunneling enhancement in nonstationary 
fields recently discovered in condensed matter physics [lO, luJ, ll^ , and can manifest itself in 
different under-barrier processes. 

For the sake of simplicity, we consider nuclear alpha decay. According to the theory 
of Gamov, the probability of the transition of alpha particle through the nuclear Coulomb 
barrier is mainly ruled by the exponential factor jisl . Q| 

W r~.exp[-A^{E)], (1) 

where 



is the classical action measured in units of h [1^. M is the mass of alpha particle, Z is the 
charge of the daughter nucleus, Rq is the nuclear radius, and the classical exit point Ra is 
determined by zero of the square root. Let us study how this probability changes when the 
decayed nuclei are placed in the beam of protons with the energy less than the barrier height 
2Ze^ I Rn. This situation rersembles the process of quantum tunneling of particles controlled 
by a wea. a.d v.yi,, ti.e e.ect_tic field ccs.de.ed in H ^ Q y H Q. 
where the specific tunneling enhancement mechanisms have been studied. The difference 
between our case and these processes is mainly in the nature of the external electromagnetic 
field (the beam of protons, in our case). The low energy projectile protons can be treated as 
a source of a pulsed electromagnetic fields interacting with the alpha decayed nuclear target. 

According to the results of papers jlfll . Illl Q|, two different regimes of tunneling are 
possible if the proton energy and its time of the under-barrier motion satisfy the necessary 
conditions. The first regime is the assistance of tunneling, when the alpha particle gains a 
part of proton energy, which can be called the positive assistance. The second one, which 
occurs when the alpha particle transfers a part of its energy to the proton, is called the 



negative assistance of tunneling. Under conditions of the negative assistance, the alpha 
particle tunnels at lower energy where the barrier is less transparent. Nevertheless, con- 
trary to any expectation, the regime of negative assistance of tunneling is unusual, since, 
under certain conditions, tunneling probability does not become exponentially small even 
for a barrier which is normally not very transparent. This phenomenon is called Euclidean 
resonance. Both mechanisms of positive and negative assistance are connected with the 
coherent multiquanta interference in the under-barrier motion. We know that the enhance- 
ment of tunneling occurs when a singularity of the nonstationary field coincides in position 
with a singularity of the classical Newtonian trajectory of the particle on the complex time 
plane. To study these processes in nuclei, we apply the semiclassical approach (based on the 
method of trajectories in complex time) developed for tunneling in nonstationary fields. 

Consider the assistance of alpha decay by the Coulomb field of incident protons when 
its energy is less than the height of the Coulomb barrier. In the absence of a proton, one 
can calculate the decay probability using formulas and (j2]). A low energy incoming 
proton interacts with the nucleus only electromagnetically, since it is stopped by the nuclear 
Coulomb field at a distance much larger than the radius of strong nuclear forces. The 
time of interactions of the proton with the decayed alpha particle (under-barrier motion 
time) is about of a characteristic nuclear time At ~ 10~^^s. The proton can excite the 
nucleus to increase the energy of emitted alpha particle E by the amount of AE during 
the interaction time At. This leads to the increase of the action in Eq. by the amoimt 
2AtAE/h (the rigorous definition of At is given below). At the same time, the energy gain 
by the excited alpha particle makes its tunneling easier due to the reduction of the action 
Aa{E) —* Aa{E + AE). Then, the resulting action for the proton induced alpha decay is 

A=^^f^ + A^{E + AE). (3) 

In condensed matter physics, this equation would describe the process of the positive photon- 
assisted tunneling with the probability exp{—A) [3]. The first term in Q results in a 
reduction of the probability due to quanta absorption and the second one describes tunneling 
in a more transparent (higher energy) part of the barrier. Since the flux of tunneling particles 
in a nonstationary field is also a nonstationary one, the expression exp(— A) relates to 
the maximal value of the tunneling probability. This maximal transition probability is 
determined by a finite value of the energy transfer AE, which provides a minimum of the 
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action A, and, hence, it is defined by the condition dA{E + AE)/dAE = 0. 

The Eq. (jS)) describes the tunnehng assisted by a given external nonstationary field. In 
our case, the role of this field is played by the proton not with a "given" but a rather with 
a "fiexible" motion affected by both the nucleus and the alpha particle. The increase of 
the value of alpha particle energy must be accompanied by the corresponding decrease of 
the proton energy value (from its initial value e down to {e — AE)). This means that the 
alpha particle and the proton participate in a cooperative tunneling motion from the nuclear 
surface to the outside of the barrier region which can be described by the joint action Aap. 
However, since the proton does not tunnel through the Coulomb barrier, the true action 



should be corrected as 
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A = A^,- ^, (4) 



where the second term (being real and negative) accounts for the phase shift of the pro- 
cess of the artificial move of the proton from outside of the barrier to the nuclear surface. 
The classical imaginary action of the proton dp can be determined from the corresponding 
Hamilton- Jacobi equation. At the limit of a weak alpha-proton Coulomb coupling, one can 
estimate 2iap/h ~ Ap{e), where Ap{e) is the proton analog of the action (j21) 



Me)-^L dr,l^-e. (5) 



Here, m is the proton mass and Zq is the charge of decayed nucleus. The action A^p is defined 
on the joint alpha-proton classical trajectory in imaginary time when both particles, with 
initial energies E (alpha particle) and e (proton), meet on the nuclear surface. The energy 
exchange between the alpha particle and the proton occurs fast and weakly contributes to 
the joint action Aap which can be written as A^p — Aa{E + AE) + Ap{e — AE). Then, the 
total action takes the form 

A = Aa{E + AE)+ Ap{e- AE)- Ap{e). (6) 

The classical time tq of the under-barrier motion, which is proportional to the derivative of 
the action with respect to energy, has the same value both for the alpha particle and for the 
proton since they move together 

2to ^ dAa{E + AE) ^ dApje - AE) 

h dAE dAE ■ ^ ^ 
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Therefore, Eq. ((Tj) determines the certain energy transfer AE which provides an extreme of 
the action ©• At the hmit of a very small energy transfer [AE -C e), the proton motion 
becomes non- "flexible" . This corresponds to the action of Eq. Q, with At = tq, which can 
be obtained from Eq. (jH)) after expansion on /\E. In this case, the tunneling motion of the 
alpha particle is affected by the nonstationary field 

V.nt{R,ir)= (8) 
\R — r[iT)\ 

where r(zr) approximately describes the classical trajectory of the free proton. The total 
action © can be written in the explicit form 



A = 



Rn nZne^\/2rn 



Then, the relation between the optimum energy transfer AE and the energy of the incident 
proton e is given by Eq. 0: 

E + AE ~ \AM) ' ^ ^ 

where we disregard the small difference between charges of the parent and daughter nuclei. It 

should be noted that generally the energy transfer AE is determined by Newtonian equations 

for the alpha particle and the proton in imaginary time, which are coupled by the interaction 

(jSI). As a consequence, the value of AE depends on the angle between directions of radial 

motions of these two particles (we consider both particles to have zero angular momenta). 

For example, AE is positive for = 180° and negative for = 0° (parallel motion). This 

means that the condition of the optimum energy transfer is fulfilled for a certain angle 

between directions of classical motion of two particles. 

One can see that the energy transfer AE during the tunneling process can be either 

positive (positive assistance of tunneling) or negative (negative assistance of tunneling). 

The latter mentioned above, is unusual, since the action can tend to zero and 

tunneling probability does not become exponentially small (Euclidean resonance) even for 

a barrier which is normally hardly transparent. Indeed, by substituting expression into 

Eq. Q we obtain 



, 27rZoe^ 2M 
A= ' 



h \ E 



f m 



1/3- 



3/2 o , ^ry -1 



If AE is negative, the proton energy e can be chosen small (see Eq. ()10j)) and the last term 
in Eq. (fTT|) may reduce the action A down to a zero value. It should be noted that the 



above equation is correct if exp{—A) <^ 1. When A becomes of the order of unity or less, 
one should use a generic formalism with the multi-instanton approach, which leads to the 
similar estimate of the action A. 

Let us give an example of the calculation of the energy transfer AE using the method 
of complex trajectories. Consider a classical parallel motion of the alpha particle and the 
proton (the angle between particles = 0) when only x-components are involved and are 
determined by the Newton equations 

M - = I • m - = - (^^\ 

dr^ Rl ^ {r, - ' dr^ (r, - R,y ^ ^ 

In the vicinity of the complex time Tq, when particles meet each other, the solutions of these 
equations have the form 

Rs Ts \ To J ^ ^ 

where Rs and are some constants to be defined. The energy AE, gained by the alpha 
particle, 

A£ = 2e=r-i^jfi (14) 

Jo [Rx - r^Y dr 

diverges close to the tq and should be regularized by the condition Rs{l — t/tqY^^ > Rq. (It 
should be noted, that in contrast to the large contribution of the diverged interaction Vint 
to AE, its contribution to the action A is not divergent and even rather small.) Then, 

\AE\ = — (^-iy\ (15) 

Rq \Rs J 



where the ratio Rs/ts satisfies the relation 



2m V r 



, RsV 2 
1 -] + 



-R.o \ ^ 1 / Rf!^ ^ 



Zo \r, 



(16) 



Considering the uraniun alpha decay as an example 



gfU + p-- IfTh + a + p (17) 

with the energy of alpha particle E = 4.678 Mev, we can fix the ratio M/m = 4, and obtain 
the parameter Rs/vg ~ 0.715. These lead to the energy transfer \ AE\ ~ 1.89 MeV and non- 
physical (negative) value of e. This means that in the case of the classical trajectory with 
= 0, the energy transfer is larger than the optimum value (see Eq. (|10|) ). The optimum 
AE (which leads to an extreme value of the action A ) corresponds to a finite of the 
classical trajectory and can be found numerically using the above scheme. 



The analysis of the expression for the action pip shows that ai e = Eq, where 

^o = E[-^) ^1.85MeV, (18) 

the energy transfer AE = 0, when the angle ^ 30°. Under those conditions, the action 
(lll|) coincides with the action of the conventional alpha decay resulting in the tunneling 
probability 

W ~ exp (4.678 MeV)] ~ e"^"-^^ - 10"^^ (19) 

The probability (fTU|) . being normalized by the nuclear attempt frequency 10^^ s^^, describes 
experimental data reasonably well [isl Il4|. 

At £ < Eq, the optimum energy transfer AE becomes finite and negative, the optimum 
angle decreases, and the action (fTT|) reduces in comparison with Aa{E). Upon reduction 
of E, the action (jllj) turns to zero at a certain proton energy er, which relates to Euclidean 
resonance. For the reaction (fTTj) = 0.25 MeV, the accompanied energy transfer is AE = 
— 1.15 MeV, and the angle between the incident proton and the emitted alpha particle is 
(j) ~ 11°. In other words, when the energy of the incident proton is about E = 0.25 MeV, the 
energy of the emitted alpha particle becomes E — \AE\ = 3.53 MeV (instead of E = 4.678 
MeV), and the energy of outgoing proton becomes £+ 1 Ai^l = 1.40 MeV. The cross-section of 
the reaction ()17j) is not exponentially small ai e = er and has a very sharp peak in the vicinity 
of the proton energy. The typical behavior of the tunneling probability as a function of the 
proton energy e is shown in Fig. ^ Positive assistance of tunneling {AE > 0) corresponds 
to the domain e > 1.85 MeV and Euclidean resonance at e = 0.25 MeV occurs at the region 
of negative assistance e < 1.85 MeV. At e > Er, the shape of the peak is proportional to 
exp[-(£: - Er)/Ae], where Ae =i lO'^MeV. 

It should be noted that at e = Er, the exit point of alpha particle is 7.3i?o and of the 
outgoing the proton is IO-Rq- The incident proton is stopped by the nuclear Coulomb field at 
the distance of about BARq, as shown in Fig.|21 Therefore, the trajectories in imaginary time 
are only a convenient language to describe the effect. In real (physical) time, the proton 
does not approach the nucleus and interacts with it solely via the Coulomb field. 

Within the exponential accuracy, the tunneling probability can be calculated by solving 
either the static problem or the equivalent dynamical one with the nonstationary field (jH)). 
As it is known, the normal tunneling through a static barrier can be described simply by 
a classical trajectory in imaginary time, which connects two classically allowed regions. 
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FIG. 1: The right side of the energy dependence (smooth enhancement) of the tunnehng probabihty 
corresponds to the domain e > 1.85 MeV and is related to positive assistance of tunnehng. The 
Euchdean resonance occurs at e = 0.25 MeV at the region of negative assistance tunnehng e < 1.85 
MeV. 

This mechanism has always been employed in the study of tunneling jiol . 20j 21, 22|. In 
the language of nonstationary field, the normal tunneling corresponds to the case when 
the tunneling particle does not absorb quanta of the nonstationary field In addition 
to this process, there is an enhanced tunneling through a barrier which corresponds, in 
the language of nonstationary field, to the absorption or emission of quanta (positive or 



negative assistance) 



12| . The description of the enhanced tunneling, in contrast to the 



normal one, does not reduce only to the classical trajectory and a narrow bundle of paths 
in its vicinity but also requires one to take into account the delicate interference of various 
paths. This can be seen from Eq. (jH), where the second part of the right side of the expression 
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FIG. 2: (a) The trajectory of the alpha particle; (b) The trajectory of the proton. These classical 
trajectories in imaginary time are only a convenient way to describe the effect. In real time the 
proton does not approach the nucleus and interacts with it solely via the Coulomb field. 

is not originated form a classical trajectory. In other words, the enhanced tunneling is a 
coherent mult i- dimensional tunneling. It should be noted that we use the method of complex 
trajectories for calculations of under-barrier processes instead of numerical solution of the 
Schrodinger equation mostly because of an insufficiency of the existing algorithms to solve 
this type of Schrodinger equation in a reasonable amount of time (see, also ref. 0]). 

The idea of stimulation of nuclear tunneling processes by incident charged particles can 
be applied in the similar way to nuclear fission processes on the basis of models with nuclear 
fragments tunneling under the action of the external varying Coulomb field. However, the 
validity of this approach for fusion reactions is less obvious and requires a further study. 



9 



In summary, protons, approaching alpha decaying nuclei, create the nonstationary 
Coulomb field, acting on the tunneling alpha particle. Due to these interactions the Eu- 
clidean resonance can appear at low proton energy and the Coulomb barrier becomes prac- 
tically transparent for the passage of the alpha particle. For example, normally, 92^!! emits 
alpha particle with the energy of 4.678 MeV. When the energy of the incident proton is close 
to its resonant value of 0.25 MeV, the energy of outgoing protons becomes 1.40 MeV and 
the energy of emitted alpha particles becomes 3.53 MeV. 
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